PASJ: Publ. Astron. Soc. Japan , 1-??, (publication date) 
© 2012. Astronomical Society of Japan. 



The Statistical Analyses of the White-Light Flares: 
Two Main Results About Flare Behaviours 
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Abstract 

We present two main results, based on the models and the statistical analyses of 1672 U-band flares. We also 
discuss the behaviours of the white-light flares. In addition, the parameters of the flares detected from two years 
of observations on CR Dra are presented. By comparing with the flare parameters obtained from other UV Ceti 
type stars, we examine the behaviour of optical flare processes along the spectral types. Moreover, we aimed, 
using large white-light flare data,to analyse the flare time-scales in respect to some results obtained from the X-ray 
observations. Using the SPSS V17.0 and the GraphPad Prism V5.02 software, the flares detected from CR Dra 
were modelled with the OPEA function and analysed with t-Test method to compare similar flare events in other 
stars. In addition, using some regression calculations in order to derive the best histograms, the time-scales of the 
white-light flares were analysed. Firstly, CR Dra flares have revealed that the white-light flares behave in a similar 
way as their counterparts observed in X-rays. As seen in X-ray observations, the electron density seems to be a 
dominant parameter in white-light flare process, too. Secondly, the distributions of the flare time-scales demonstrate 
that the number of observed flares gets a maximum value in some particular ratios, which are 0.5 or its multiples, and 
especially positive integers. The thermal processes might be dominant for these white-light flares, while non-thermal 
processes might be dominant in the others. To reach better results for the behaviour of the white-light flare process 
along the spectral types, much more stars in a wide spectral range, from the spectral type dK5e to dM6e, must be 
observed in the white-light flare patrols. 

Key words: methods: data analysis — methods: statistical — stars: activity — stars: flare — stars: individ- 
ual^ Dra) 



1. Introduction 

The flare process has not been perfectly understood yet, 
although several studies on this subject have been carried 
out since the first flare was detected on the Sun by R. C. 
Carrington and R. Hodgson on September 1, 1859 (Carrington 
1859; Hodgson 1859). The flare activity of dMe stars is mod- 
elled based on the processes of Solar Flare Event (Gershberg 
2005; Benz & Giidel 2010). To understand flare event of dMe 
stars, the levels of the flare energy have been looked over 
in many studies. Gershberg & Shakhovskaya (1983) demon- 
strated that there are some differences between energy levels 
of different stars and star groups. They demonstrated that the 
flare energy levels of the Orion association stars are located at 
the highest levels. The Pleiades stars are located just below 
the Orion association stars. Finally, the stars in the Galactic 
field are located below all of them. The levels seem to be ar- 
ranged according to the different ages. However, the differ- 
ences among the different stars and star groups might be due 
to the saturation level of white-light flares detected from UV 
Ceti type stars. The X-ray and radio observations indicate that 
some parameters of magnetic activity can reach the saturation 
(Gershberg 2005; Skumanich & McGregor 1986; Vilhu et al. 
1986; Doyle 1996; Garcfa-Alvarez et al. 2008; Jeffries et al. 
2011). Dal & Evren (2011a) recently revealed some clues for 
the saturation in the white-light flares. They modelled the dis- 
tributions of flare-equivalent durations versus flare total dura- 



tion. In the models, it is seen that flare-equivalent durations 
can not be higher than a specific value and it is no matter how 
long flare total duration is. According to Dal & Evren (201 la), 
this level defined as the Plateau parameter in the models is 
an indicator of the saturation level for the white-light flares in 
some respects. 

To understand flare event of dMe stars, the flare time-scales 
are as important as emitting flare energy. There are some pa- 
rameters (such as B, n e and geometry of flaring loop) to deter- 
mine the border of the time-scales in the flare events. These pa- 
rameters affect the time-scales of heating and cooling of the re- 
gions, where the flare events occur (Temmer et al. 2001 ; Reeves 
& Warren 2002; Imanishi et al. 2003; Pandey & Singh 2008). 
According to Imanishi et al. (2003), the flare decay time (r^) 
is firmly correlated with B and the electron density (n e ), while 
the flare rise time (i>) is proportional to a larger magnetic loop 
length (£) and smaller B values. Reeves & Warren (2002) con- 
cluded that the observed light curve is modelled with using 
multi-loop instead of a single-loop. However, some discrepan- 
cies were observed in the decay phases and especially toward 
the end of the flare light curves. In contrast to that assumed by 
Imanishi et al. (2003), Pandey & Singh (2008) found that there 
is a sustained heating during the decay phases of the flares. 
Moreover, Favata et al. (2005) demonstrated that short flares 
occur in a single loop, while longer flares are generally arcade 
flares and some heating is driving their decay phases. In addi- 
tion, there are some relations between the flare rise and decay 
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times. Using the event asymmetry index (A ev ), Temmer et al. 
(2001) statistically analysed solar H a flares, and found some 
relations between flare decay and rise time-scales. Kay et al. 
(2003) demonstrated that there is a linear correlation between 
flare rise and decay times. 

In the first part of this study, the flares of CR Dra were mod- 
elled, and some parameters were derived from the model. We 
compared CR Dra with other UV Ceti type stars in respect to 
these parameters. In this point, CD Dra takes an important 
place, because it fills a gap in the B-V range. CR Dra is a 
flare star (Gershberg et al. 1999), which is classified as a dou- 
ble or multiple star in the SIMBAD database. The orbital pe- 
riod is 4.04 yr and the semimajor axis of the orbit is 0".148. 
The total mass of the system is reported as a value between 
1.0 Mq and 1.8 M . The distance of the system is given be- 
tween 17.4 pc and 20.7 pc (Tamazian et al. 2008). CR Dra 
is a metal-rich stars and a member of the old disk population 
in the galaxy (Stauffer & Hartmann 1986; Veeder 1974). The 
photometric magnitude difference between the components of 
the system was measured as l m .8 in V-band by Tamazian et 
al. (2008). The flare activity of CR Dra was discovered for the 
first time by Petit (1957). Apart from the flare activity, there are 
some studies for the sinusoidal-like variations at out-of-flares, 
such as Mahmoud (1991);Mahmoud (1993); Anderson (1979). 
However, it is not clear that the star exhibits any rotational 
modulation. Mahmoud (1991) found that CR Dra exhibits a 
sinusoidal-like variation with period of 16 d in V-band, while 
the period found from B-band is 5 d . The author mentioned 
that the amplitude of the variation found in B-band was lower 
that the observational standard deviation. On the other hand, 
Mahmoud (1993) found that the period of the sinusoidal-like 
variation at out-of-flare is 8 d .776 in V-band. 

In the second part of this study, a large flare data containing 
1672 U-band flares, which were obtained in this study and col- 
lected from 42 different studies in the literature, were analysed 
to test whether any relation exists between the flare rise and 
decay times for the white-light flares, or not. 

2. Observations and Data 

2.1. Observations 

Observations were acquired with the High-Speed Three- 
Channel Photometer (HSTCP) attached to the 48 cm 
Cassegrain-type telescope at Ege University Observatory. 
Observations were performed in two different manners. We 
used a tracking star in the second channel of the photometer 
and carried out the flare observations only in standard Johnson 
U-band with exposure times between 2 and 10 seconds and 
time resolution of 0.01 seconds. The second type of observa- 
tions was used for determining whether there was any variation 
out-of-flare. We also observed CR Dra once or twice a night, 
when the star was close to the celestial meridian. These ob- 
servations were made with the exposure time of 10 seconds 
for each band of standard Johnson BVR system. Considering 
the technical properties of the HSTCP given by Meistas (2002) 
and following the procedures outlined by Kirkup & Frenkel 
(2006), the mean average of the standard deviations of obser- 
vation times was computed as 0.08 seconds for the U-band ob- 
servations. It was found to be 0.20 seconds for the multi-band 



observations. The same comparison stars were used for all ob- 
servations. The basic parameters of all the program stars (such 
as standard V magnitudes and B-V colours) are given in Table 
1. Although the program and comparison stars are so close 
in the plane of the sky, differential extinction corrections were 
applied. The extinction coefficients were obtained from ob- 
servations of the comparison stars on each night. Moreover, 
the comparison stars were observed with the standard stars in 
their vicinity and the making use of differential magnitudes, 
in the sense of variable minus comparison, were transformed 
to the standard system using the procedures outlined by Hardie 
(1962). The standard stars are listed in the catalogue of Landolt 
(1983). Heliocentric corrections were applied to the obser- 
vation times. The mean averages of the standard deviations 
are m .015, m .009, m .007 and m .007 for the observations 
acquired in standard Johnson UBVR bands, respectively. To 
compute the standard deviations of observations, we use the 
standard deviations of the reduced differential magnitudes in 
the sense comparisons (CI) minus check (C2) stars for each 
night. There is no variation in the standard brightness compar- 
ison stars. The differential magnitudes in the sense of compar- 
ison minus check stars were carefully checked for each night. 
The comparison and check stars were found to be constant in 
brightness during the period of observations. 

In general, 20 U-band flares were detected in observations 
of CR Dra. The star was observed 5\38 in 2005 and 29' 1 .10 
in 2007. All the flares were detected in U-band observations of 
2007. Apart from flares, no variation was found in BVR bands 
out-of-flare activity. Some samples of detected flares are shown 
in Figure 1 . Several parameters, such as the flare rise time, de- 
cay time, amplitude and equivalent duration, were computed 
from the light curve of each flare. The procedure used in the 
calculations is similar to the methods found in the literature 
(Gershberg 1972; Moffett 1974). The procedure schema of cal- 
culations are shown in Figure 2 for two flares, which are seen 
in both top and bottom panels of Figure 1. In calculations, we 
firstly separated each flare light curve into three parts. One of 
them is the part indicating the quiescent level of the brightness 
before the first flare on each nigh. The brightness level without 
any variations (such as a flare or any oscillation) was taken as a 
quiescent level of the brightness of this star. To determine this 
level, we used the standard deviation of each observation point, 
considering the mean average of all the observation points un- 
til this last point. If the standard deviations of the following 
points get over the 3a level, this point was taken as the begin- 
ning of a flare. The quiescent levels of each star were deter- 
mined from all the observation points before the first flare on 
each night. It must be noted that, although a different quies- 
cent level was found for each observing night, their levels are 
almost the same considering the 3a levels of observing con- 
dition. As seen in Figure 2, we fitted this level with a linear 
function (/i(x) and gi(x)), and then, using this linear func- 
tion, we computed the flare equivalent duration, flare ampli- 
tude and all the flare time-scales (rise and decay times). The 
part of the light curve above the quiescent level was also sepa- 
rated into two sub-parts. First of them is the impulsive phase, 
in which the flare increases. Second one is the decay phase. 
The impulsive and decay phases were separated according to 
the maximum brightness observed in this part. It must be noted 
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that some flares have a few peaks. In this case, the point of the 
first-highest peak was assumed as the flare maximum. To de- 
termine the flare time-scales, we fitted the impulsive and decay 
phases with the polynomial functions. In Figure 2, f%(x) and 
52 0*0 are lne fi ls of the impulsive phases for the given sam- 
ples, while fs(x) and 53 (x) are the fits of decay phases. The 
best polynomial functions were chosen according to the cor- 
relation coefficients (r 2 ) of fits. To determine the beginning 
and end of each flare, we computed the intersection points of 
the polynomial fits with the linear fit of the quiescent level and 
their standard deviations. In this study, the intersection points 
were taken as the beginning and end of each flare. The flare 
rise time (r r ) was taken the duration between the beginning 
and the flare maximum point. For example, the duration from 
the first intersection between fx (x) and ji (x) to the second one 
between /2(x) and fz(x) was taken as the flare rise time (r r ) 
for this flare. In the same way, the flare decay time (t<;) was 
taken the duration between the flare maximum point and the 
flare end. In Sample 1, the duration from the second intersec- 
tion between /2(x) and /a(x) to the third one between fz{x) 
and /i(x) was taken as the flare decay time (r^). The height 
of the observed-maximum points from the /1 (x) linear fit was 
taken as the amplitudes of this flare. The same procedure was 
used for each flare, and the Maple 12 (Monagan et al. 2008) 
software was performed in all calculations. 

Apart from the flare time-scales and amplitudes, using 
Equations (1) and (2) taken from Gershberg (1972), the equiv- 
alent durations and energies of all the flares were computed: 

P = J[{Iflare-Io)/h]dt (1) 

where Iq is the flux of the star in the observing band while 
in the quiet state, and Ifi are is the intensity at the moment of 
flare. For instance, the parameter P is equal to about total area 
between the /i(x), /2(x) and fs(x) for Sample 1. 

E = PxL (2) 

where E is the energy, P is the flare-equivalent duration in the 
observing band, and L is the intensity in the observing band 
while the star is in the quiet state. 

For each observed flare, the HJD of flare maximum mo- 
ment, flare rise (r r ) and decay time (r^) (s), flare total durations 
(s), ratio of Td/r r , flare-equivalent duration (s), flare amplitude 
(mag), and their energies (erg) were calculated. Those param- 
eters are given in Table 2. 

2.2. Data Used In the Analyses 

The number of flares used in the analyses are 1672 from 
which 20 flares are detected in the observations of CR Dra in 
this study. 534 flares were obtained from other stars of this 
project. The remaining 1118 flares were collected from 41 dif- 
ferent studies in the literature. Thus, the data used in the anal- 
yses were contained the parameters of 15 different UV Ceti 
stars. In Table 3, the distribution of the data taken from the 
literature is listed versus references ordered by year. 

In the analyses of flare activity, especially in the OPEA mod- 
els, 554 flares obtained in this project were only used. This is 
because the analyses need to be comprised of parameters de- 
rived with the same method and from the flares detected with 
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the same optical system. Otherwise, some artificial variations 
and differences can occur between the data sets. To avoid this 
problem, instead of all 1672 flares, just 554 flares detected 
in this project, which comprised parameters derived with the 
same method and the same optical system, were used. On the 
other hand, all 1672 flares were used in the analyses of the 
distribution of the flare numbers versus the ratio of flare de- 
cay time to rise time. Determining and computing of the time- 
scales of the rise and decay phases are almost the same in all 
the studies. 

3. Analyses 

3.1. Flare Activity and the One-Phase Exponential 
Association Models 

The distributions of the equivalent durations in the logarith- 
mic scale versus flare total durations were derived in order to 
test whether there are any upper limits for the distributions of 
the equivalent durations (logP u ) of the flares detected in the 
observations of CR Dra. Although the flare energy is generally 
used to analyse in the literature, the flare equivalent durations 
were used to analyse in this study. This is because of the lu- 
minosity parameter (L). As seen from Equations (2), the flare 
energy (E) depends on the luminosity parameter (L) with flare 
equivalent duration (P). However, the luminosity (L) is differ- 
ent for each star. Although there are small differences among 
the masses of M dwarfs, the luminosities of the two M dwarfs, 
whose masses are very close to each other, can be dramati- 
cally different due to their position in the H-R diagram. This 
means that the computed energies of flares are very different 
from each other, even if the light variations of the flares occur- 
ring on these two stars are absolutely the same. Because of this, 
the equivalent durations (P) were used instead of energy (E) 
in this study. The equivalent duration parameter (P) depends 
only on flare power. 

In order to model this distribution, first of all, the best curve 
was estimated with SPSS V17.0 software (Green et al. 1999) 
and GraphPad Prism V5.02 software (Motulsky 2007). The 
regression calculations showed that the Exponential function 
is the best model. Besides, the tests done with GraphPad 
Prism V5.02 software revealed that especially the One Phase 
Exponential Association function (hereafter OPEA) (Motulsky 
2007; Spanier & Oldham 1987) given by Equation (3) is the 
best exponential function among all others. The OPEA model 
of the distributions of the equivalent durations in the logarith- 
mic scale versus flare total durations was derived with using the 
Least-Squares Method with GraphPad Prism V5.02 software. 

y = y + (Plateau - 50) X (1 - e~ k x x ) (3) 

There are some important parameters derived from this func- 
tion, which have some clues about the flaring loop and the con- 
dition inside the loop. One of them is yo, which is the lower 
limit of equivalent durations for observed flares in the loga- 
rithmic scale. In contrast to yo, the parameter of Plateau is 
the upper limit. The value of yo depends on the quality of ob- 
servations as well as flare power, while the value of Plateau 
depends only on power of flares. In fact, according to Equation 
(2), the value of Plateau depends only on the energy of flares 
occurring on the star. The parameter k in Equation (3) is a con- 
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stant depending on the x value, which is the total flare time. 
Apart from them, one of the most important parameter is the 
Half — Life parameter. This one is half of the first x val- 
ues, where the model reaches the plateau values for a star. The 
Half — Life parameter was also computed from the derived 
model. The derived OPEA model is shown in top panel of 
Figure 3, and the parameters of the model are listed in Table 4. 
In the middle panel of Figure 3, the model derived for CR Dar 
is compared with the models derived for other six UV Ceti type 
stars (Dal & Evren 201 la; Dal & Evren 201 lb). In the bottom 
panel of the figure, the Plateau values are compared for these 
stars. 

Considering the flares in the Plateau phases of the OPEA 
model, the mean average of maximum flare-equivalent du- 
rations was computed to test whether the Plateau value is 
statistically acceptable, or not. The Independent Samples t- 
Test (hereafter t-Test) (Wall & Jenkins 2003; Dawson & Trapp 
2004) was used for the calculation. Thus, the value of Plateau 
was tested using another statistical method. Although the mean 
value computed by the t-Test is expected to be close to the 
Plateau value of the OPEA model, it is clear that there can be 
some difference between these two values. This is because the 
OPEA model depends on all the distributions of the x values 
from the beginning to end, while the mean value computed by 
the t-Test depends only on equivalent durations of flares in the 
Plateau phase. The results obtained from the t-Test analyses 
are also listed in Table 4. 

Examining the flares detected in U-band observations of CR 
Dra, it was computed that the maximum flare rise time obtained 
from these 20 flares is 1967 s, while the maximum flare total 
duration is 4955 s. 

3.2. Distribution of the Flare Numbers versus the Ratio of 
Flare Decay Time to Flare Rise Time 

Several parameters were computed from the light curves of 
all these flares detected in this project. Then, using some statis- 
tical methods, the relations among themselves of all the param- 
eters were analysed with both SPSS V17.0 (Green et al. 1999) 
and GraphPad Prism V5.02 software (Motulsky 2007). During 
the analyses of the distribution of flare amplitudes versus the 
ratio of flare decay time to rise time, a remarkable gathering 
was seen for some ratios of flare decay time to rise time, which 
values are several times of the values 0.5 and especially the 
positive integer numbers, such as 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 
3.5, 4.0, 4.5, etc. These accumulations are seen in Figure 4. 
As seen from the figure, the ratios of flare decay times to rise 
times are especially positive integers for many flares with dif- 
ferent amplitudes. 

To examine this unusual unexpected case, first of all, the 
ratio of flare decay time to rise time were computed for each 
flare. The values of / r r are listed for CR Dra flare in Table 2. 
Considering the errors listed in Table 2, the mean error of / r r 
values was found to be ±0.042 for the flares detected from CR 
Dra. However, in this project, 554 U-band flares are detected 
in total. Taking all 554 flares into account, the mean error of 
T~d/T r values was found to be ±0.026. As also seen from Table 
2, although there are some flares, whose Td/r r errors are dra- 
matically larger than the computed mean errors, the total num- 
ber of these flares is small enough to enough to be neglected. 



Apart from the flare detected in this project, we also collected 
1118 flares from the literature and combined them with 554 
flares. The sources of 1 1 18 flares are listed in Table 3. It must 
be noted that unfortunately we could not compute the errors of 
Td/Tr values for the flares taken from the literature due to ab- 
sence of errors of r r and Td values in the sources. Secondly, 
the best histograms were determined using SPSS V17.0 and 
GraphPad Prism V5.02 software for these data. Table 5 shows 
the parameters use to produce the histograms. The analyses in- 
dicate that the best statistically acceptable interval ratio length 
should be 0.05 in the interval ratio length. 

The derived numbers of the flares in intervals of 0.05 were 
plotted versus the ratio of flare decay time to rise time. Then, 
the histogram was derived for 554 flares detected in this 
project. The ratio distribution is shown in Figure 5, while the 
results of the analysis is listed in Table 5. As seen from the 
figure, the ratio is equal to 1.0 in 0.05 ratio length for 73 flares, 
and it is 2.0 for 39 flares and also 3.0 for 29 flares. In brief, 
the ratios of flare decay time to rise time are especially equal 
to 1.0 or its multiples as the positive integers for 204 flares 
among 554 flares. Besides, the ratio is 0.5 for 18 flares. In to- 
tal, the ratios of flare decay times to rise times are equal to 0.5 
or its multiples, such as 1.5, 2.5, 3.5, 4.5, etc. The incidence 
of the flares, whose ratio of their decay times to rise times is 
1.0 or positive integers, is 36.82% over 554 flares detected in 
this project. It is 15.70% for the flares, whose ratios are 0.5 or 
its multiples, such as 1.5, 2.5, 3.5, 4.5, etc. Considering that 
the positive integers are acceptable as the multiples of 0.5, the 
incidence of the flares, whose ratios of their decay times to rise 
times are 0.5 or its multiples, such as 1.0, 1.5,2.0,2.5,3.0,3.5, 
4.0, 4.5, etc., is 52.53% over 554 flares. 

In order to test whether this remarkable result obtained from 
the data of this project might be occurred due to any selec- 
tion effects, almost all the U-band flare data given in the liter- 
ature were collected from 41 different studies. Thus, 1118 U- 
band flares were collected apart from 554 U-band flares of this 
project. In fact, there are more 170 flares given in literature, but 
the decay or/and rise times were not given. Because of this, the 
ratios of flare decay times to rise times were not computed for 
these ones, and they could not be used. Following the method 
mentioned above, the same examinations were done for these 
1118 U-band flares. The result is seen in Figure 6, while the 
results of the analysis is listed in Table 5. According to the sta- 
tistical analyses, the incidence of the flares, whose ratio of their 
decay times to rise times is 1.0 or positive integers, is 16.37% 
over 1118 flares collected from the literature. It is 8.23% for 
the flares, whose ratios are 0.5 or its multiples, such as 1 .5, 2.5, 
3.5, 4.5, etc. Considering that the positive integers are accept- 
able as the multiples of 0.5, the incidence of the flares, whose 
ratios of their decay times to rise times are 0.5 or its multiples, 
such as 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, etc., is 24.60% over 
these 1118 flares. 

Finally, both the data obtained in this project and the data 
collected from the literature were combined together in order to 
reach rather better result. In brief, a data set containing 1672 U- 
band flares in total was obtained. Following the same method, 
the distribution was derived and statistically examined. The re- 
sult of the examinations is shown in Figure 7, while the results 
of the analysis is listed in Table 5. Consequently, considering 
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1583 U-band flares, the incidence of the flares, whose ratio of 
their decay times to rise times is 1 .0 or positive integers, was 
found as 23.09%. It was computed as 10.83% for the flares, 
whose ratios are 0.5 or its multiples, such as 1.5, 2.5, 3.5, 4.5, 
etc. Considering that the positive integers are acceptable as the 
multiples of 0.5, the incidence of the flares, whose ratios of 
their decay times to rise times are 0.5 or its multiples, such as 

1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, etc., is 33.91% over these 
1583 flares. 

4. Results and Discussion 

4.1. Flare Activity and the One-Phase Exponential 
Association Models 

The distributions of flare-equivalent durations versus flare 
total duration were modelled by the OPEA function expressed 
by Equation (3) for CR Dra. As it is seen from Figure 3, this 
function demonstrates that the equivalent durations of the flares 
detected in the observations of CR Dra can not be higher than a 
specific value regardless of the length of the flare total duration. 
In order to test whether the Plateau parameter derived from 
the model is statistically acceptable, using the t-Test, the mean 
average of equivalent durations was computed for the flares lo- 
cated in the Plateau phase in the OPEA models. The mean 
average of equivalent durations was found to be close to the 
Plateau values derived from the OPEA models. The Plateau 
parameter derived from the OPEA model was identified as an 
indicator of the saturation level for the white-light flares by Dal 
&Evren (2011a). 

According to Dal & Evren (2011a), two stars, EV Lac and 
EQ Peg, have almost different level of the Plateau value. 
However, the observations of CR Dra have demonstrated that 
the variation of the Plateau value is different from that shown 
by Dal & Evren (2011a). In fact, the shape of the Plateau 
value variation is remarkably different from the conclusions of 
Dal & Evren (2011a). EV Lac and EQ Peg are not different 
from the others. CR Dra flares indicate that the variation of 
this parameter has a trend, and both EV Lac and EQ Peg are on 
the trend. The observations of CR Dra also demonstrated that 
white-light flare activity could be affected by some parameters 
of the flare processes as well as the flares detected in X-ray 
or radio observations. The figure taken from Katsova et al. 
(1987) demonstrated that electron density varies along B-V in- 
dexes (Gershberg 2005), and its variation trend is in contrast 
shape with that of the Plateau values. The analyses of flares 
detected in X-ray have revealed some important points of the 
processes. Considering the parameters of the Solar Flare Event, 
the suspected parameters could be one or several of va, B, R 
or n e . The suspected parameters must be B or/and n e in the 
atmosphere of the stars during the flare event. According to 
Katsova et al. (1987) it could be n e . On the other hand, Doyle 
(1996) suggested that it can be related to some radiative losses 
in the chromosphere instead of any saturation in the process. 
However, Grinin (1983) demonstrated the effects of radiative 
losses in the chromosphere on the white-light photometry of 
the flares. According to Grinin (1983), the negative H opac- 
ity in the chromosphere causes the radiative losses, and these 
are seen as pre-flare dips in the light curves of the white-light 
flares. 



The other parameters derived from the OPEA model are the 
flare time-scales. The maximum flare duration was computed 
as 4955 s, while the maximum flare rise time was computed as 
1967 s from the detected flares of CR Dra. The Half — Life 
parameter was found to be 191.40 s from CR Dra flares. All 
these time-scales were found to be in agreement with those 
found for other stars by Dal & Evren (2011a); Dal & Evren 
(2011b). It is well known that the decay time-scales of the 
X-ray light curves is related to the length of the loop (Haisch 
1983; van den Oord et al. 1988; van den Oord & Barstow 
1988; Reale et al. 1997; Reale & Micela 1998; Reale et al. 
1988). According to the results found by Favata et al. (2005), 
short flare durations could be an indicator of single loop geom- 
etry, while long flare durations could reveal arcade flares. If 
the same case is valid for the white-light flares observed from 
UV Ceti type stars, the lengths of the flaring loops on the sur- 
face of CR Dra are generally larger than those on AD Leo, EV 
Lac, EQ Peg according to the time-scales obtained from CR 
Dra. However, the general sizes of the flaring loops are almost 
equal to loops occurring on VI 005 Ori. In the general perspec- 
tive, the decreasing of the observed maximum time-scales (the 
flare rise time and the flare total duration) of flares reveals that 
the lengths of the flaring loop are decreasing toward the later 
spectral types among UV Ceti type stars. On the other hand, 
the flare mechanism of CR Dra easily reaches the saturation in 
shorter time than the stars such as AD Leo, EV Lac, EQ Peg. 
In this point, CR Dra and VI 005 Ori have a same nature. 

It must be noted that it is a debated issue whether the white- 
light flares can be affected from the parameters such as ua, 
B, R and n e , which generally affect the coronal loops and the 
flare emissions in the corona (Benz & Giidel 2010; Hawley & 
Fisher 1992; Abbett et al. 1999; Hawley et al. 2003; Alfred et 
al. 2006). However, the results found in this study reveal that 
there are some similarities between the white-light flares and 
their counterparts observed in the X-ray or radio bands. 

Consequently, the observations of CR Dra have demon- 
strated that the Plateau values vary in a trend from earlier 
spectral types to the later. However, the trend is not in the 
shape previously offered by Dal & Evren (201 la). Although it 
is a debated issue, the shape of the trend in the Plateau vari- 
ation versus B-V indexes is in agreement with the trend in the 
distribution of n e versus B-V indexes. However, much more 
stars between the spectral type dK5e and dM6e should be ob- 
served in order to reach more reliable result. Especially, more 
stars around CR Dra should be observed. 

4.2. Distribution of the Flares Numbers versus the Ratio of 
Flare Decay Time to Flare Rise Time 

The most remarkable result was found in the analyses of the 
distribution flare amplitudes versus the ratios of decay times to 
rise times. As it can be seen from Figure 4, although the error 
bars of some flares are a bit large, the flares are gathering in the 
specific ratios. According to the distribution seen in the fig- 
ure, many flares seem to prefer the value of 0.5 or its multiples, 
such as 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, etc., in the ratios 
of decay times to rise times. To test whether these accumula- 
tions in these ratios are statistically meaningful, the histograms 
were derived and analysed to find statistically the best one. The 
derived best histograms are shown in Figures 5, 6 and 7. As 
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seen from the figures, the statistical analyses indicated that the 
best width of bars is 0.05 for histograms. The mean error is 
found to be ±0.042 for t^/tv of the flares detected from CR 
Dra. Although the error is close to the estimated width value 
of bars, the mean error is lower. In addition, when the mean 
error computed over all 554 flares is ±0.026. In this case, the 
histograms are statistically acceptable. 

The histograms demonstrated some definite results. First of 
all, most of the flares prefer the value of 1.0 or its multiples 
as positive integers. The incidence of them over 1672 flares is 
23.09%. The incidence of the flares, which prefer the value of 
0.5 or its multiples, such as 1.5, 2.5, 3.5, 4.5, etc., is 10.83%. 
The positive integers can be taken as the multiples of the value 
0.5. In this case, the incidence of the flares, whose ratios of 
the decay time to rise time are 0.5 and its multiples, reach the 
value of 33.91% over 1672 U-band flares. However, as it is 
stated in the Section 3.2, the incidences are found a bit differ- 
ent from the analysis of each data set. This must be because 
of some small differences in the methods used to compute the 
flare time-scales, such as decay or rise times. In this project, 
all the parameters were computed with considering the quies- 
cent level of the brightness of a star. These levels were com- 
puted from the part of the observations without any flares or 
any variations for each night separately. On the other hand, in 
the literature, the first points of the flare beginnings were taken 
as a quiescent level of the brightness of observed star in some 
studies. Unfortunately, the first point of the flares does not al- 
ways indicate the quiescent level of the brightness. Taking the 
first points of the flare beginnings as a quiescent level comes 
some difference to computed decay and rise times. Using the 
data collected from different studies might cause a bit different 
incidences. 

In the case of many flares, preferring the ratios as 0.5 
and its multiples might be related with some parameters in 
the flare event process, which affect the flare time-scales 
such as the flare rise or decay times. The flare decay time 
(Td) is firmly correlated with B and n e , while the flare rise 
time (tv) is proportional to a larger I and smaller B values 
(Temmer et al. 2001; Reeves & Warren 2002; Imanishi et al. 
2003; Pandey & Singh 2008). The magnetic loop length (£) 
also depends on plasma electron density (n e ) and plasma tem- 
perature (Yokoyama & Shibata 1998; Shibata & Yokoyama 
1999; Shibata & Yokoyama 2002; Yamamoto et al. 2002). As 
it is seen, the preferred ratios of the decay times to rise times 
is mainly related with the electron density (n e ) in the flaring 
loop and quite a bit with magnetic field strength {B) of the 
loop. In fact, with some assumptions, Imanishi et al. (2003) 
has shown the relations between flare time-scales (flare decay 
and rise times) and plasma electron density (n e ) and the recon- 
nection factor (Ma) with Equations (A5) and (A8) given by 
them. 

It is seen from Figures 5, 6 and 7, there are two type flares. 
One of them is the flares discussed above, which are shown 
by filled and open circles in the figures. Other type flares are 
shown by small points in the figures, which are neither the mul- 
tiples of 1.0 nor the multiples of 0.5. These two types must 
be the slow and fast flare types (Gurzadian 1988). Gurzadian 
(1988) indicated that thermal processes are dominant in the 
processes of slow flares, which are 95% of all flares observed 



in UV Ceti type stars. Non-thermal processes are dominant in 
the processes of fast flares. In fact, the r r and Td parameters of 
the fast flares can take random values according to themselves 
(Gurzadian 1988; Dal & Evren 2010). Thus, the small points, 
whose ratios are neither the multiples of 1.0 nor the multiples 
of 0.5, must indicate the fast flares. 

Consequently, the distribution of the number of flares versus 
the ratio demonstrates that (1) the number of observed flares 
gets maximum in the ratio, which are 0.5 or its multiples, and 
especially positive integers. As a result of this, Td/r r is usu- 
ally equal to 0.5 or its multiples, and mostly a positive integer. 
(2) The flare numbers are dramatically decreasing toward the 
larger values of the ratio. (3) Some flares do not prefer specific 
values of the ratio, thus their ratios of the decay times to rise 
times are neither the multiples of 1.0 nor the multiples of 0.5. 
(4) The thermal processes might be dominant in the processes 
of the flares preferring specific ratios, while non-thermal pro- 
cesses might be dominant in the others. 
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Fig. 1. Three samples for the observed flare light curves in U-band. The flare shown in top panel was observed in April 1 8, 2007, while flares shown in 
middle and bottom panels were detected in May 13, 2007. In the figures, filled circle represents observations, while dashed line represents the quiescent 
level of the brightness determined for each night. 
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Fig. 2. Two samples for computing of the flare parameters. In both figures, the filled-black circles represent observations in the quiescent level of the 
brightness. The blue-squares represent the impulsive phase of the flares, while the red-diamonds represent the main (decay) phases of the flares. The 
black lines (fi(x) and gi(x)) represent the linear fits of the observations in the quiescent level, while the blue lines (f2(x) and g2(x)) represent the 
linear fits of the observations in the impulsive phase. The red lines (f:i(x) and gz (x)) represent the polynomial fits of the decay phases of the flares. 
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Fig. 3. Comparison of CD Dra's The OPEA with other stars. Top panel: The OPEA mode] derived from the distributions of flare-equivalent duration on 
a logarithmic scale vs. flare total duration for CR Dra. Filled circles represent equivalent durations computed from flares detected from CR Dra. The 
line represents the model identified with Equation (3) computed using the least-squares method. The dotted lines represent 95% confidence intervals for 
the model. Middle panel: The OPEA model (solid line) of CR Dra is compared with the models (dotted line) derived from observations on other 6 stars 
reported by Dal & Evren (201 la); Dal & Evren (201 lb). Bottom panel: The Plateau parameters vs. the B-V index of the CR Dra and other 6 stars taken 
from Dal & Evren (2011a); Dal & Evren (201 lb) are shown. In the panel, open symbols represent the parameters taken from Dal & Evren (201 la); Dal 
& Evren (201 lb), while filled circle represents the parameter of CR Dra. 
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Fig. 4. A sample for the distribution of the flare amplitudes versus the ratio of flare decay time to flare rise time. Instead of any amplitude variation in a 
regular trend, there is a remarkable gathering in some particular ratios of flare decay time to rise time. 
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Fig. 5. The distribution of the number of the flares derived in intervals of 0.05 ratio length versus the ratio of flare decay time to rise time. Top panel: the 
distribution is shown for ratios between 0.00 and 58.35. All the flares seen in this figure were obtained in this project. The value of 58.35 is the largest 
ratio obtained from the flares detected in this work. Bottom panel: the distribution is shown for the interval between 0.00 and 10.00 for clarity. In both 
panels, the filled circles represent the flares, whose ratios of flare decay times to rise times are equal to 1 .0 or its multiples as the positive integers. The 
open circles represent the flares, whose ratio is 0.5 or its multiples, such as 1.5, 2.5, 3.5, 4.5, etc. The small points represent the flares with other ratios, 
which are neither the multiples of 1.0 nor the multiples of 0.5. 
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Fig. 6. The distribution of the number of the flares derived in intervals of 0.05 ratio length versus the ratio of flare decay time to rise time. Top panel: 
the distribution is shown for all ratios from 0.00 to 206.40. All the flares seen in this figure were collected from the literature. The value of 206.40 is the 
largest ratio obtained from these flares. Bottom panel: the distribution is shown for the interval between 0.00 and 10.00 for clarity. In both panels, all the 
symbols used in this figure are the same with the symbols used in Figure 4. 
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Fig. 7. The distribution of the number of the flares derived in intervals of 0.05 ratio length versus the ratio of flare decay time to rise time. Top panel: the 
distribution is shown for all ratios from 0.00 to 206.40. The data plotted in the figure are a combination of data from this work and those taken from the 
literature. The value of 2006.40 is seen to be the largest ratio obtained from all the flares. Bottom panel: the distribution is shown for the interval between 
0.00 and 10.00 for clarity. In both panels, all the symbols used in this figure are the same with the symbols used in Figure 4. 
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Table 1. Basic parameters for the targets studied, comparison (CI) and check (C2) stars. 



Program Stars 


V (mag) 


B-V (mag) 


CR Dra 


9.997 


1.370 


CI: HD 238545 


10.195 


0.582 


C2: HD 238552 


10.123 


1.021 



Table 2. All the parameters computed for each flare detected in flare patrols of CR Dra in this study. All the flares were obtained in the U-band. 



HJD of Flare 


Rise 


Decay 


Total 




Equivalent 


Flare 


Flare 


Maximum 


Time 


Time 


Duration 


Td/T r 


Duration 


Amplitude 


Energy 


(+24 00000) 


(s) 


(s) 


(s) 




(s) 


(mag) 


(erg) 


54208.47638 


60±0.5 


96±0.6 


156±0.8 


1.60±0.03 


20.3403±0.0007 


0.188±0.003 


- 9.5898 E+30 


54208.48038 


110±1.1 


2090±2.1 


2200±2.4 


19.00±0.2 


183.0304±0.0024 


0.518±0.008 


- 2.4743 E+32 


54209.53320 


305±1.0 


394±2.3 


699±2.5 


1.29±0.06 


75 .2661 ±0.0006 


0.270±0.004 


- 3.5485 E+31 


54210.53322 


204±0.5 


341±1.1 


545±1.2 


1.67±0.01 


42.23 18±0.0005 


0.255±0.004 


- 1.9911 E+31 


54210.54189 


48±0.6 


288±1.5 


336±1.6 


6.00±0.10 


33.3097±0.0012 


0.299±0.004 


- 1.5704 E+31 


54211.32701 


48±0.6 


33±0.4 


81 ±0.7 


0.69±0.02 


11.6477±0.0004 


0.329±0.005 


-5.4915 E+30 


54211.48255 


120±0.4 


1245 ±0.6 


1365±0.8 


10.38±0.04 


112.4808±0.0009 


0.496±0.007 


- 1.0080 E+32 


54227.37541 


15±0.2 


75±0.9 


90±0.9 


5.00±0.12 


14.5103±0.0005 


0.484±0.007 


-6.8411 E+30 


54227.49538 


15±0.2 


15±0.2 


30±0.3 


1.00±0.02 


5.0416±0.0002 


0.422±0.006 


-2.6598 E+30 


54227.50538 


1044±1.5 


1260±2.1 


2304±2.6 


1.21±0.01 


133.2189±0.0013 


0.258±0.004 


-8.1931 E+31 


54234.38692 


364±0.4 


741 ±0.4 


1106±0.5 


2.03±0.01 


27.8203±0.0010 


0.277±0.004 


- 1.3116E+31 


54234.50690 


373±0.6 


787±0.9 


1160±1.1 


2.11±0.01 


164.2169±0.0027 


0.488±0.007 


-7.7422 E+31 


54235.47382 


1967±1.1 


2988±1.3 


4955±1.7 


1.52±0.01 


145.4049±0.0031 


0.503±0.008 


- 1.3288 E+32 


54276.37291 


45 ±0.5 


15±0.2 


60±0.6 


0.33±0.01 


8.5775±0.0003 


0.291 ±0.004 


~ 4.0440 E+30 


54276.39774 


90±1.1 


210±2.5 


300±2.7 


2.33±0.06 


51.0005±0.0018 


0.342±0.005 


-2.4045 E+31 


54276.40208 


15±0.2 


60±0.7 


75±0.7 


4.00±0.10 


10.8760±0.0004 


0.230±0.003 


-5.1277 E+30 


54276.40399 


105±0.3 


225 ±0.3 


330±0.4 


2.14±0.01 


50.4296±0.0018 


0.247±0.004 


- 2.3776 E+31 


54276.41007 


120±0.2 


255 ±0.4 


375±0.4 


2.13±0.01 


65.9255±0.0013 


0.256±0.004 


-3. 1082 E+31 


54276.41927 


30±0.2 


30±0.2 


60±0.2 


1.00±0.01 


12.2865±0.0004 


0.465 ±0.007 


- 5.7927 E+30 


54276.43185 


30±0.2 


15±0.2 


45±0.2 


0.50±0.01 


6.6471 ±0.0002 


0.441 ±0.007 


- 3.1339 E+30 
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Table 3. Summary of the reference list and the data used in the analyses. 



[Vol. , 



Reference 


Observed 


Spectral 


Number 


Reference 


Observed 


Spectral 


Number 


of Data 


Star 


Tvne* 


of Flare 


of Data 


Star 


Tvne* 


of Flare 


MacConnell 11968") 

IVXtlt/ V^VJllll\_.ll ( x y \J 1J I 


AD 1 eo 


M4.5Ve 


21 


Panov # Tsvetkov (4 981) 

1 tlllVJ V OC 1 o V \^LIVVJ V 1 1 yUl 1 


EV Lac 


M4.5V 


11 


MacConnell CI 968) 


VI 054 Onh 


M3.5Ve 


1 1 


Panov et al 11982a) 

X HI III V V- L Ll 1 ■ I 1 y U ^ CI / 


EV Lac 


M4.5V 


7 


Herr & Brcich (1969) 


DO Cep 


M4.0V 


10 


Panov etal. (1982b) 


AD Leo 


M4.5Ve 


22 


Cristaldi & Rodono 11970) 


AD Leo 


M4.5Ve 


2 


Osawa et al 11973") 

v.y j tl VV ll v l Ll 1 ■ 1 ly / J 1 


AD Leo 


M4.5Ve 


14 


Cristaldi & Rodono f 1 Q71 ^ 

^ 1 1 old 1111 L-X_ XVVJV1VJ11VJ \ ± ^/ 1 ± J 


EV Lac 


M4.5V 


3 


Panov etal (1983") 

A CI 1 III V \^ L £11. 1 1/U J 1 


EV Lac 


M4.5V 


12 


Cristaldi & Rodono 11972) 


AD Leo 


M4.5Ve 


1 


Herr & Frank 11983") 

1 1 V_ 1 1 Lv X 1 111 1 lv \ 1. ^/ KJ *J } 


AD Leo 

/ \ 1 / V^V 1 


M4.5Ve 


6 


Cristaldi & Rodono 11973a) 


EV Lac 


M4.5V 


3 


Pettersen et al (1 984) 

X V- LLL1 . > V- 1 1 V- L Ll 1 . \ 1 .y U T 1 


AD Leo 

/ \ L^* 1 / V-V 1 


M4.5Ve 


82 


Rodono (T973) 

xvvjixvjiivj I 1 y / J 1 


EV Lac 


M4.5V 


4 


Panov et al 11985") 

X tlllVJ V \^ L til. 1 l/UJ J 


EV Lac 


M4.5V 


20 


Kanoor et al 11971") 

XVtlLJVJWl \^L til. \ L J 1 ^3 } 


AD 1 eo 

Aly 1 f \^\ } 


M4.5Ve 


2 


Pptfprsen et al 11 986) 

X V UL1 .)V 11 \^L til. 1 1 y yj\J } 


AD I eo 

i \ i_y 1 - A^v } 


M4.5Ve 


99 


Cristaldi & Rodono C1973W 


EV Lac 


M4.5V 


1 


Tsvetkov et al 11 986a) 

1 □ VL L1VW V Ll 1 . 1 l/UUily 


EV Lac 


M4.5V 


5 


Moffett 11974) 

1V1\J11 L'LL \ ly / r 1 


Wolf 424 


M5 


1 1 


Tsvetkov et al 11 986b) 

1 □ L1VW V Ll 1 . 1 1 y UUU / 


EV Lac 


M4.5V 


17 


Moffett 11974") 

1V1U11 L LL \ 1 ^/ 1 r / 


UV Cet 


M6.0V 


52 


Herr & Onie 11987") 

i iv. 1 1 L_X- v.y iy i v_. \ \_ j \j t } 


AD Leo 


M4.5Ve 


14 


Moffett 11974) 


EV Lac 


M4.5V 


21 


Herr & Canuto 11988) 

1 L V_ 1 1 V_\- V.^ Ll ll Ll LU \ X ^/ U U 1 


DO Cen 

l^* v.y v.^ V- 1 ) 


M4.0V 


5 


Moffett (1974) 


EO Pes 


M3.5 


42 


Herr & Charache (1988) 


AD Leo 


M4.5Ve 


6 


Moffett 11974)) 

1V1U11 V- 1 . 1 11/ / r 1 1 


YZ CMi 


M4.5V 


51 


Panov etal 11988") 

x hi ivy V til . \ X ^/ U U } 


EV Lac 


M4.5V 


22 


Moffett 11974") 

1VXVJ11C-LL ^17 ' '/ 


YY Gem 


dMle 


14 


Pettersen # Snndland 11991) 

1 \_- 1 1 V_ 1 .)V 11 LX- O L11HJ1 till V.I I X .y J X 1 


V577 Mon 

V *J 1 1 1.VXVJ11 


M4.5V 


72 


Moffett 11974) 

1VXV_J11\_- LL 1 x y 1 ~ J 


AD Leo 


M4.5Ve 


8 


Antovetal 11991") 

nn lu v ll tii . i i y y i j 


AD Leo 

ni/ i 


M4.5Ve 


4 


Moffett (1974) 


CN Leo 


M6.0V 


102 


van den Oord et al. (1996) 


YZ CMi 


M4.5V 


8 


Osawa et al. (1974a) 


YZ CMi 


M4.5V 


7 


Letoetal. (1997) 


EV Lac 


M4.5V 


193 


Kapoor&Sinvhal(1974) 


AD Leo 


M4.5Ve 


4 


Tovmassian et al. (1997) 


EV Lac 


M4.5V 


20 


Osawa et al. (1974b) 


AD Leo 


M4.5Ve 


4 


Panov et al. (2000) 


YZ CMi 


M4.5V 


7 


Moffett (1975) 


AD Leo 


M4.5Ve 


27 


Dal & Evren (2010) 


EV Lac 


M4.5V 


98 


Sanwal (1975) 


YZ CMi 


M4.5V 


1 


Dal & Evren (2010) 


AD Leo 


M4.5Ve 


110 


Nicastro (1975) 


DO Cep 


M4.0V 


22 


Dal & Evren (2010) 


V1054Oph 


M3.5Ve 


40 


Ichimura & Shimizu (1978) 


YZ CMi 


M4.5V 


7 


Dal & Evren (2010) 


EQ Peg 


M3.5 


73 


Contadakis et al. (1979) 


EV Lac 


M4.5V 


23 


Dal & Evren (2011a) 


V1005Ori 


MOVe 


41 


Contadakis et al. (1980) 


EV Lac 


M4.5V 


13 


Dal & Evren (2011b) 


VI 285 Aql 


M3.0V 


83 


Avgoloupis et al. (1980) 


EV Lac 


M4.5V 


5 


Dal (2011) 


DO Cep 


M4.0V 


89 



* All the spectral types were taken from the SIMBAD database. 



Table 4. The parameters derived from the OPEA model and the results of the t-Test analyses. 



Parameter Value ±Error 

Parameters of The OPEA Function: 
Plateau = 2.146±0.065 
2/o = 0.673±0.088 
k= 0.003621 ±0.000584 

Span= 1.473±0.097 
Half-Life= 191.40 
r 2 of the model = 0.94 

The t-Test Analysis: 
Mean= 2.163±0.037 
Std. Deviation = 0.082 



No. 



Usage of pas j 00 . els 
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Table 5. The values of the best histograms, which were estimated by SPSS V17.0 and Prism V5.02 software. 



Parameters 


In This Project 


Literature 


Combined 




(Figure 5) 


(Figure 6) 


(Figure 7) 


Total Number of Flares = 


554 


1118 


1672 


Minimum Value of tj)t t = 


0.07 


0.09 


0.07 


25% Percentile = 


1.00 


2.02 


1.51 


Median Value of t^/tv = 


2.00 


4.79 


3.50 


75% Percentile = 


4.31 


11.13 


8.81 


Maximum Value of tJt t = 


58.33 


206.39 


206.39 


Mean Value of t c iIt t = 


4.16 


9.35 


7.63 


Std. Deviation = 


6.03 


14.01 


12.21 


Std. Error = 


0.26 


0.42 


0.30 


Lower 95% CI of mean = 


3.66 


8.53 


7.05 


Upper 95% CI of mean = 


4.66 


10.18 


8.22 



